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A solution has recently been given [Jones and Gallet, 1962] to the problem of representing 
the complex variations of ionospheric characteristics on a worldwide scale, including their 
diurnal variation, by numerical analysis of ionospheric data as measured at the stations, 
without prior hand operations. Whereas the paper referred to above deals with ''how to 
make numerical maps," the present paper is concerned with ''how to use a numerical map" 
once it has been obtained. Included are (1) a precise definition of a numerical map, (2) 
efficient computing procedures for applying numerical maps, including a general method for 
computing worldwide (or polar) contour maps in either universal time or local mean time, 
and (3) a large selection of graphical illustrations computed automatically from a numerical 
map. 



1. Introduction 

In a recent paper by Jones and Gallet [1962] a 
solution was given to the problem of representing 
the complex properties of ionospheric characteristics 
on a worldwide scale, including their diurnal varia- 
tion, b}' numerical analysis of ionospheric data as 
measured at the stations, without prior hand opera- 
tion.^ The solution referred to above is general 
enough to be applied to any ionospheric characteris- 
tic, and, in fact, has already been applied to the 
following: the critical frequency (/o^2), the 3000-kin 
maximum usable frequency factor (F2-M3OOO), the 
maximum electron density (A^max), the height of 
Amax (Hmax), ^ud tlic quartcr thickness of a layer 
(SCAT) [see Wright, Wescott, and Brown, 1961]. 
With only shght modification, the methods could 
also include variations with height above the surface 
of the earth. 

Ionospheric representations (or maps) obtained 
from this solution are given in a numerical form and 
are therefore referred to as numerical maps (sec. 2). 
Whereas the first reference [Jones and Gallet, 1962] 
deals with ''how to make numerical maps/' the pres- 
ent paper is concerned with ''how to use a numerical 
map'' once it has been obtained. 

A discussion of methods for applying numerical 
maps is given in section 3. Included are efl&cient 
computing procedures for evaluating a numerical 
map and a general method for computing world- 
contour maps in either universal time or local mean 
time. A wide selection of such graphical representa- 
tions is given in section 4. Graphical displays of 
these types are needed for solving problems of HF 
radio propagation by manual operations and, more- 
over, serve a useful purpose in describing the complex 
properties of the ionosphere. One should recognize, 



1 For a brief summary of the problem and solution, the reader can refer to 
[Jones and Gallet, 1960]. 



however, that when electronic computing equipment 
is available, it is more accurate and usually more 
efficient to use numerical maps directly within the 
computer than to read and interpolate from contour 
maps. 

As a means of illustration, the numerical map used 
in the present paper is for the monthly median of 
/o7^2, September 1958. The types of illustrations 
given here could also be given for any other charac- 
teristic, and for predicted numerical maps for future 
months, since these maps have exactly the same 
mathematical form. We begin in the following 
section by giving a precise definition of a numerical 
nuip. 

2. Numerical Maps r(X,^,t) 

The term, numerical map, is used to denote a 
function T(\,6,t) of three variables, latitude (X), 
longitude (6), and time (^), whicli represents an iono- 
spheric characteristic, including its diurnal variation, 
on a worldwide scale. In order to apply such a 
function once it has been obtained, one must know 
its exact form and the range of definition of each 
independent variable. We begin with the latter: 



\= geographic latitude: 



-90°^X^ 90°, (1) 



e= geographic longitude: 0° ^ (9 ^ 360°, (2) 

(6 in degrees east of Greenwich) 

t=^local mean hour angle: — 180° ^ t ^ 180°, (3) 
^=(15/^—180°), where A=local mean time 
(LMT) in hours. 

For example, at noon LMT, h=l2 and t=0^. 
Universal time (UT) can therefore be introduced by 
means of the relation 



T=t-d, -180°^ T^ 180°, 



(4) 
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where T is the universal hour angle. 

The general form of T(X,d,t) is a 'Tourier'' time 
series 



T(\,9,t)=--ao(\,d)+^[aj{\,d) cos jt + bj(\, 6) sin jt] 

(5) 

where H denotes the number of harmonics retained 
for the diurnal variation and where the ^Tourier^' 
coefficients aj(\,d) and bj{\,d) vary with the geo- 
graphic coordinates.^ Such variations are repre- 
sented in the form 



k=0 

k=0 



3=0,1,..., H, 



H, 



(6) 



where the geographic functions Gk(\,d) are given in 
table 1. As can be seen the functions 0^(^,6) consist 
of three main groups, each having different types of 
longitudinal terms present. Thus a specific set of 
the functions GjcO^,d) is specified by assigning the 
values qo, ^i, and ^2 which are, respectively, the 
highest powers of sin X for terms involving (1) no 
longitudinal variation, (2) first order longitudinal 
variation, and (3) second-order longitudinal varia- 
tion. Equivalently, one could also specify the values 
ko, kij and Z:2=^ taken by the index k at the end of 
each of the three groups of the G;,{\,d), respectively. 
These values are related to the q^ values by 

ko=qo and ^-^=^^-1+2(^^+1) for r=-l,2. (7) 

A numerical map T{\,d,t?) can therefore be written 
more explicitly in the form 

T(\,d,t)=^Do,jcG,{\,d) 

k=0 

+i: -(l j:D,j,,G,x\e)'] cos jt 

,G,{'K,e)'jsinjt\ 



L/v-o 



(8) 



which is defined by a relatively small table of 
coefficients Dsk- An example of such a set of coeffi- 
cients is given in table 2 for the characteristic 
/o^2 monthly median for September 1958. As can 
be seen from table 2, the number of harmonics 
retained is H=S and the parameters required to 
specify the functions Gjc(\,6) are given by ko=ll, 
ki=^37, and k2=K=4:d. Hence from (7) we have 
for the highest powers of sin X in each of the three 
groups of the Gk(\,d), qo=n, ^i=12, and ^2 = 3. 
The blanks in table 2 indicate zero coefficients. 



2 The systematic and well-defined geographic variations of these Fourier co- 
efTicients are illustrated by an atlas of graphs for monthly median ^^2 for four 
seasonal months and for high and low years of solar activity [Jones, 1962J. 





Table 1. 


Geographic functions GkO^,d) 




Main latitudinal 




Mixed latitudinal and longitudinal variation 


variation 


First order in longitude 


Second order in longitude 


k 


Gk(\,d) 


k 


Gk(\,e) 


k 


Gk(\,d) 



1 
2 


1 

sin X 
sin2 X 


ko+l 
/co+2 

/co+3 
ko+4. 


cos X cos d 
CCS X sin d 

sin X cos X cos 
sin X cos X sin 9 


fci+l 
ki-\-2 

/ci+3 
/ci+4 


cos2 X COS 29 
cos2 X sin 29 

sin X cos2 X cos 20 


k(^ 


sin«oX 


sin X cos2 X sin 26 














kx-1 
ki 


sin^i XccsX cos 9 
sin?! X cos X sin 


K-1 
K 


sin 92 X cos2 X cos 20 
sin«2 X cos2 X sin 20 



Before considering methods for applying numerical! 
maps in section 3, we insert here a brief word ofj 
caution. The functions G)c(K,d) (in table 1) are] 
not orthonormal relative to the coordinates (Xt,^i)| 
of the ionospheric stations; hence one cannot obtain 
a best representation (in the sense of least squares) : 
of lower degree by merely truncating the series! 
in (6).^ When a numerical map is originally | 
produced, aj(\,d) and bj{\,6) are represented in I 
terms of orthonormal functions, corresponding to 
the set of stations available for the particular month. 
Such series are then truncated at optimum places 
determined by the noise and are subsequently 
expressed in the form of (6) for convenience and 
simplicity in practical applications [Jones and 
Gallet, 1962, ch. 2]. 

3. Methods for Applying Numerical Maps 

For applying numerical maps, a number of useful 
methods and techniques have been developed which 
eliminate slow and tedious hand work as much 
as possible and which make efficient use of automatic 
computing and plotting equipment. The discussion 
of such methods given in the present section is 
not intended to be exhaustive in nature but rather 
a starting place for developing further methods for 
applying numerical maps. Included are such topics 
as (1) efficient methods for evaluating r(X,^,^), (2) 
a general procedure for computing contour maps, 
and (3) a numerical check on the coefficients Dsk- 
The methods described below have been programed 
for use on several large-scale digital computers, the 
IBM 704 and 7090 and the CDC 1604. 

3.1. Evaluation of r(X, 6, t) 

Using the coefficients in table 2, one can compute 
the value of/o7^2 monthly median for September 1958 
at any location on the earth and for any desired 
instant of time. This is the basic operation in most 
applications of numerical maps. The evaluation of 
r(X, 6, t) for fixed values of X, 6, and t can be performed 
in the straightforward manner indicated by (8), 



3 the Gk(\,0) in table 1 are used, since they are the simplest functions which are 
linear combinations of the classical surface spherical harmonics [Byerly, 1893]. 
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TABLE 2 






















TIME VARIATION 














Harmonic 





1 


2 


3 


4 






^^ 





1 


2 


3 4 


5 


6 


7 


8 







1.2566249E 01 


2,26259G9E 00 


1.2521171E 00 


-1.2278355E 00 1.3986979E 00 


-3.7681434E-01 


-6.8212361E-01 


2.4211725E-01 


-3.7045519E-01 








1 


2.3950864E-01 


-1.2906564E 00 


-2.5089329E-01 


1.6727505E GO 1.0199610E 00 


1.7951671E-02 


-1.3118478E 00 


-7.2732570E-02 


1.1056117E-01 










2 


-1.1465794E 01 


5.7717085E-01 


6.1885204E 00 


3.0063888E 00 -4,6062881E 00 


3.5320769E 00 


-4.4741024E 00 


3.0824802E-01 


3.9941964E 00 










3 


1.8361369E 01 


4.1578991E 01 


-6,2927405E 00 


-2.7043140E 00 -3.2191761E 00 


-3.2875664E 00 


2.2687291E 01 


9.5334392E-02 


3.6307619E-01 










4 


4,6198351E 01 


-1.5204600E 01 


8.6436892E-01 


-1.3568193E 01 8.1038432E 00 


-1.9840365E 01 


3.9411273E 01 


-6.1881873E-01 


-1.5809304E 01 








I 


5 


-1.2754594E 02 


-2.7716534E 02 


3.5377018E 01 


3,6834899E-01 -8.3601906E 00 


2.7900419E 01 


-1.2042979E 02 




-2.2943625E 00 








6 


-1.3831402E 02 


1.4315429E 01 


-2.5665253E 01 


3.9922174E 01 -9.0052624E 00 


3.7436317E Ol 


-1.0610687E 02 




2.3262825E 01 










7 


3.0226488E 02 


6.9731399E 02 


-7,6536765E 01 


-1.6047472E 01 3.3546848E 01 


-7.0617263E 01 


2.7192087E 02 




1.3728741E 00 










8 


1.6182682E 02 


2.0255935E 01 


3.1032109E 01 


-4.8789722E 01 8.2568223E 00 


-2.81085>41E 01 


1.1904225E 02 




-1.1227530E 01 










9 


-3.0370904E 02 


-7,4465693E 02 


6.9262616E 01 


3.6984787E Ol -3.7459259E 01 


7.1432502E 01 


-2.7368133E 02 














10 


-6.3623919E 01 


-2.2118568E 01 


-1.2634591E 01 


2.0873692E 01 -4,3315212E 00 


7.3103423E 00 


-4.7419598E 01 










2 

o 




11 


1.1015161E 02 


2.8443719E 02 


-2.2136436E 01 


-2.0425142E 01 1.4921944E 01 


-2.53846i7E 01 


1.0112482E 02 










12 


-1.7257335E-01 


5.1166203E-02 


1.0165302E-01 


-7.0306255E-02 -2 . 91 73936E-0 1 


-6.6573759E-03 


2.0101654E-01 










13 


3.1468502E-01 


-1.4744951E-01 


-1.2851838E-01 


8.6758951E-02 -2 .6830614E-0 1 


-8.9503505E-02 


-1.4615241E-01 










J— 




14 


2.1917868E 00 


-2.6453310E 00 


2.0596705E 00 


-7.2709031E-01 2,2399525E 00 


-5.5057557E-03 


-3.6947747E-01 










< 




15 


-1.1110747E 01 


-1.1914294E 01 


5.0979509E 00 


7.6853372E-0i -2.6008389E 00 


-1.7334515E 00 


-1.0453669E 00 










tr 




16 


-1.9544060E 01 


-1.2653360E 01 


3.8076664E 00 


3.0899957E 00 5.5101019E 00 


-1.3538277E 00 


-2.8346492E 00 










< 




17 


1.3633266E 01 


8.7443910E 00 


-1.0040889E 01 


-8.0549182E-01 ^.4819739£ 00 


2.2136988E 00 


2.2643434E 00 










> 




18 


-1.0187599E 01 


2.6621571E 01 


-2,7641896E 01 


4.7570975E 00 -4,2233717E 01 


-1.32l434^c 00 


1.7484305E 00 










—1 




19 


2.3785661E 02 


1.8263532E 02 


-1.2897653E 02 


-3.5782082E 00 7,4267178E Ol 


1.2384460E 01 


2.6761956E 00 










< 




20 


2.1337638E 02 


1.4022231E 02 


-2.6366782E 01 


-3.2841640E 01 -3.6472735E Ol 


3.6292006E 00 


7.3703185E 00 










o 




21 


-2.3088949E 02 


-5.3263182E 01 


1.1863364E 02 


1.3797593E 00 -4,9261112E Ol 


-5.6528990E 00 


-5.4238232E 00 










X 




22 


-2.9804478E 01 


-1.2770266E 02 


1.4492764E 02 


-9,9443520E 00 2.1907150E 02 


3.9350583E 00 


-1.5457780E 00 










< 




23 


-1.2235357E 03 


-7,9008117E 02 


6.9895783E 02 


-4.5919808E 00 -4.5075892E 02 


-2.5599535E Ol 


-1.4800732E 00 










II 


24 


-8.2631350E 02 


-5.4231029E 02 


-1.2023700E 01 


1.5074860E 02 1.1935836E 02 


-2.3942029E 00 


-5.1717052E 00 










O 


25 


1.1058288E 03 


2.4813321E 01 


-4.1681761E 02 


-7.015170^E 00 2.0773483E 02" 


3.9299649E 00 


3.3373286E 00 








ay 
en 


o 

LtJ 




26 


2.2010015E 02 


2.7810015E 02 


-3,8857270E 02 


1.2801645E 01 -4.5419638E 02 


-2.8032897E 00 












27 


2.6505530E 03 


1,3954952E 03 


-1.5475186E 03 


2.9632579E 01 1.1057194E 03 


1.5516487E 01 










H-i 






28 
29 
30 
31 
32 
33 
34 
35 
36 
37 


1.4967556E 03 

-2,3196936E 03 

-3.6275302E 02 

-2.5672058E 03 

-1.2937514E 03 

2.2341890E 03 

1.8481805E 02 

9.1761112E 02 

4,3247125E 02 

-8.0814595E 02 


1.0250898E 03 

3.1640388E 02 

-2.7396303E 02 

-1.0686661E 03 

-9.6218758E 02 

-5.7781020E 02 

I.OOIUOOE 02 

2.9208834E 02 

3.5444336E 02 

2.8306951E 02 


2.7741387E 02 

6.7919226E 02 

4.6111605E 02 

1.5359677E 03 

-4.6372998E 02 

-5.4909590E 02 

-1.9298043E 02 

-5,6471215E 02 

2,2504081E 02 

1.8094877E 02 


-3.1285267E 02 -2.2276641E 02 
2.7396283E 01 -4.1083528E 02 

-1.5688963E 01 4,0905625E 02 

-4,1375109E 01 -1.1920315E 03 
3.0056282E 02 2.2661473E 02 

-3.4611450E 01 3.8728193E 02 
6.0569722E 00 -1.3366369E 02 
1.8240356E 01 4,6979535E 02 

-1.0803830E 02 -9.3319020E Ol 
1.4428984E 01 -1.4080027E 02 














38 


4.75794285-01 




-2.4874791E-01 


1.5660529E-01 1 . 12 18772E-0 1 
















39 


-1.5646780E-01 




-2,5391683E-01 


-9.9137233E-02 -1 .00022 77E-0 1 


















40 


1.5809660E 00 




-8,6182244E-01 


1.5531334E-01 9.6352962E-01 
















III 


41 


-8.5030371E-01 




-1.7133347E-01 


2.3344462E-01 1 . 36 768 52E-0 1 
















42 


-1.1362623E 00 




7,2939227E-01 


-3,8717439E-01 


















43 


-8.5501653E-01 




3.0101221E 00 


-4.6i95528E-0l 


















44 


-2.1345755E 00 






-1.7114903E 00 


















45 


1.6638957E 00 






1.4107662E-01 














_J 
< 
O 

X 

< 

q: 


2 
O 

< 
0^ 


Harmonic 


5 


6 


7 


8 






K^ 


9 


10 


11 12 


13 


14 


15 


16 







1.5436883E-01 


2.5i34453E-0l 


-1.6932602E-01 6 . 1690475E-02 


-6.5109502E-02 


-1.2627347E-01 


9.1846^21E-02 


-4.5376377E-02 




CD 
O 
UJ 


^ 


I 


1 


-2,5472259E-02 


7.6121460E-03 


-7.7820108E-03 -1 . 5589729c-02 


2.1036162E-03 


-2.3577334E-03 


-6.4546649E-03 


4,9245555E-04 






> 




2 


-1.7009369E-01 


-2.4425167E-01 


1.8239843E-01 -8 .3470014E-02 


7.5380495E-02 


1.4439629E-01 


-1.0264139E-01 


8.6725935E-02 






I - Mam lot 


itudinal variation. Mixed latitudinal an 


d longitudinal variation II -First order in longitude, 


III - Second order in longitude 




1 






Notation- For 


gach entry the number given by the f 


rst eight digits and sign is multiplied by the power of ten defined by the last two digits and 


Sign. 1 






CO 


EFFICIENTS D^^ DEFINING 


THE FUNCTION r(X ,^,t) 


FOR MONTHLY MEDIAN 


fQF2(Mc/s) 




' 










SEPTEMBER 


958 









provided the coefficients Vsk in table 2 are given.^ 
However, it is frequently required to compute the 
value of r(X, 6^ t) (a) at many locations for a fixed 
instant of time, or (b) at a fixed location for several 
different instants of time. In such cases consider- 
able savings in computer cost can be achieved by 
using one of the following sets of procedures: 
(a) Problem: 

Given a fixed value ^o of the local mean hour 
angle. 

Compute r(X, 6, Iq) for several points (X, d). 
Procedure: 

Step 1. Evaluate cos jt^ and sin jt^ j=l, 2, 

Step 2. Compute : 

H 

Dk{to)=Do^k+^ [D2j,k cos jto+D2j-i,k sin jto] 



for each ^=^0, 1, . 
coefficients Dsk. 

Step 3. Evaluate Gj,{\ 6), k^O, 1, . 
for one of the given points (X, 6), 



, K using step 1 and the 
, K 



Step 4. Compute r(X, d, ^o)=S DMG,(\, d). 

Step 5. Repeat steps 3 and 4 for each remain- 
ing point (X, 6). 

Remark: The economy of this technique lies in 
the fact that steps 1 and 2 do not have to be 
repeated so long as ^o remains fixed. 

(b) Problem: 

Given fixed values of longitude ^o ^^nd universal 
hour angle Tq, 

Compute r (X, do, to), where to=do-{-To, for several 
values of X. 
Procedure: 

Step 1. Follow steps 1 through 5 in (a) above, 
letting ^0=^0+^0 and giving the fixed value do 
throughout. The points (X, B) in step 3 become 

(c) Problem: 

Given fixed values of latitude Xo and longitude 

Compute r(Xo, ^o, i) for several values of the 
local mean hour angle t. 
Procedure: 



Step 1. Evaluate Gk{\o, do), 



Step 2, Compute: 

K 

and 



^-0, 1, 



j=o,h 



K. 



,H, 



K 

^jO^Oj Oo)^=^ D2j-i,]cGic(>^o, ^o) , 
k=0 

i=i,2 

* As was pointed out in section 2, the parameters H=S, A:o=ll, /:i^ 
50=11, gi = 12, and 72 =3 can be determined from table 2. 



... .,H. 

--^7,k2=K=i5, 



Step 3. Evaluate sin jt for j=l, 2, 
and cos jt for ./=1, 2, 
for one of the given values of t. 
Step 4' Compute: 



, H 



H 

r(Xo, do, t)=ao(\o, ^o)+S \(^j(\, <9o) cos jt 

+ bj(Ko,6o) sm jt]. 

Step 5. Repeat steps 3 and 4 for each 

required value of t. 

Remark: The economy of this technique lies in 
the fact that steps 1 and 2 are not repeated so 
long as the point (Xo, ^o) remains fixed. 

3.2. General Method for Computing Contour Maps 

In section 4 several types of contour maps of 
monthly median /oi^2 are illustrated: (a) world and 
polar maps in universal time, (2) world maps in 
local mean time, and (3) maps in latitude and local 
mean time for fixed longitudes. Each of these maps 
has been computed using the coefficients Dsic in 
table 2 and the general procedure described here. 
The procedure is of interest not only from the view- 
point of the applications illustrated in the present 
paper, but also since it can be applied to a wide 
variety of problems in which physical or other 
variations have been mapped in a numerical form. 

The basic problem is the following: to locate 
(for a given time t) all points along certain meridian 
lines at which the function T{\,d,t) assumes any one 
of a given set of values Iq, h, . . ., Iq. For simplicity 
we assume the values Ig to be equally spaced and 
increasing and hence defined by 



/, 



/o</,. 



q=0, 1, . . ., Q. 



(9) 



A solution to this basic problem is given in the 
appendix. By repeated application of this solution 
one can compute any of the three types of contour 
maps mentioned above. The procedures are as 
follows. 

a. World Maps in Universal Time 

Let To denote a given universal hour angle (sec. 2) 
and let 



.=360«(g) 



m=0, 1, 



M, 



and 



tm — lo\dm 



m=0, 1, ..., M. 



-0, 1, 



(10) 

(in 

M we 



Then for each pair {dm,tm) m^ 
compute all values of X satisfying 

T(\,dm,tm) = I, for some ^=0, 1, ..., Q, (12) 

using the method given in the appendix. By 
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choosing M sufficiently large, the points on the map 
can be brought closely enough together to define 
clearly each contour. In practice M can be nuide 
as large as we please, but for purposes of economy 
it is desirable to make it as small as can be allowed. 
For most ionospheric characteristics, ikf== 24 is 
large enough. 

b. World Maps in Local Mean Time 

Let ^0 denote a given local mean hour angle (sec. 2) 
and let dm be given by (10). Using the method in 
the appendix, we compute (for each dm) all values of 
\ satisfying 

T{\,dm,to) = I, for some q=0, 1, ..., Q. (13) 

Again the number M is chosen as described above. 

c. Maps in Latitude and LMT for Fixed 6 

Let ^0 denote a given longitude and let 

tm=lS0''[2mJM-l], 7n = 0, 1, ..., M (14) 

denote a set of local mean hour angles. Again 
using the solution given in the appendix, we can 
locate all values of X such that 

T{\,do,tm) = I, for some (j = 0,l, ..., Q (15) 

for each value of m. As before, we must choose M 
large enough so that the pouits on the map clearly 
define each contour. M=24 is sufficient for most 
ionospheric characteristics. 

Before considering the graphical illustrations in 
section 4 computed by the above procedures, we 
describe briefly a useful numerical clieck which can 
be made on the coefficients Dsk before each applica- 
tion. 

3.3. A Numerical Check on the Coefficients IDsk 

In the normal operations of handling, sorting, and 
reproducing decks of punclied cards and converting 
to magnetic tapes and reading into a computer, 
there is always the possibility that a card will be 
misplaced or that a number will be transferred or 
reproduced incorrectly. For this reason it is recom- 
mended that certain periodic checks be made when 
applying the coefficients Dsk defining a numerical 
map. 

One such test that can be performed automatically 
by the computer is the following. Given a set of 
correct values of the function T(X,d,t) for a few values 
of the independent variables X, 6, and t, one can test 
all of coefficients Dsk by recomputing the given 
values of r(X,^,0. If the functional values agree 
(to, say, six digits) then the coefficients can be 
used with confidence. As an example we give a 
set of functional values in table 3 that can be used 
to test the coefficients in table 2. 



Table 3. Functional values of r(X,^,t) for testing coefficients 
I^sk for monthly median foF2 {Mc/s), September 1958 (see 
table 2) 

(Fixed: X = 50.0°, 0=90.0°) 



LMT 


t 


Y{\e,t) 






Mc/s 


GO 


-180° 


6. 8834498 


01 


-165° 


6. 6014013 


02 


-150° 


6. 2867521 


03 


-135° 


6. 0181661 


04 


-120° 


5. 7592264 


05 


-105° 


6. 0961668 


06 


- 90° 


7. 6478070 


07 


- 75° 


9. 6392061 


08 


- 60° 


10. 9863478 


09 


- 45° 


11. 7429891 


10 


- 30° 


12. 2211317 


11 


- 15° 


12. 3494360 


12 


0° 


12. 1784082 


13 


15° 


11. 8577095 


14 


30° 


11. 4744272 


15 


45° 


11. 1244725 


16 


60° 


10. 8307340 


17 


75° 


10. 5763241 


18 


90° 


10. 2428354 


19 


105° 


9. 6362065 


20 


120° 


8. 8870622 


21 


135° 


8. 2365329 


22 


150° 


7. 6430851 


23 


165° 


7. 1626558 



4. Graphical Representations From 
Numerical Maps 

As was mentioned above, many hiteresting and 
useful results can be prockiced from a numerical 
map. The basic application, using tlio coefficients 
Dsk, is to compute the value of T(\,d,t} — i.e., of the 
ionospheric characteristic^ — ^at any desired location 
and instant of time in the three-dimensional space 
of latitude, longitude, and local mean time.^ For 
the purpose of illustration, a schematic diagram of 
tliis space is given in figure 1. In this diagram 
several types of (two-dimensional) cuts and (one- 
dimensional) lines are shown, on each of which a 
graphical representation of the ionospheric cliarac- 
teristic can be produced from its numerical map. 
These cuts and lines, summarized in table 4, should 
be considered only a suggestive sample and not an 
exhaustive set. 

Table 4. Summary of cuts and lines in the three-dimensional 
space of X,6,t 

(F=fixed; V= variable) 





Latitude 


Longitude 


LMT 


UT 


Surfaces: 










Cut I 


V 


V 


F 




Cutn 


V 


F 


V 




Cut ni 


V 


V 




F 


Lines: 










Line A 


F 


F 


V 




LineB 


F 


V 


F 




Line C 


V 


F 


F 




LineD 


F 


V 




F 



5 Simple procedures for making this application were given in section 3.2. 
These can be performed with such speed by an electronic computer that very 
extensive applications are practical. 
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i_ _ 


B 




\ 


L__._. 



SilTl! 




Figure 1. Three-dimensional space of \,B,i. 

X = geographic latitude, ^ = geographic longitude, 
i( = local mean hour angle, 
T^ universal hour angle. 



In this section, examples are given of graphical 
representations for each cut and line in table 4, such 
representations being computed directly from the 
coefficients in table 2. For each line in the three- 
dimensional space (fig. 1), the graphical representa- 
tion is a simple two-dimensional graph computed 
from V{\d,t) in the obvious manner. For cuts, 
however, the representations take the form of con- 
tour maps on which selected values of V{\,6,t) are 
traced throughout a plane region. Such maps are 
computed by solving high degree polynomial equa- 
tions in sin X for points lying on selected meridian 
lines (sec. 3.2). Since T{\,d,t) is defined continuously 
throughout the space, points on the contours can be 
obtained as close together as desired. The following 
are illustrations of such graphs and contour maps. 

4.1. World Maps in Universal Time 

World maps of monthly median ^0^2 in universal 
time (Cut III) are illustrated in figures 2 and 3. 
Such maps are probably the most useful for problems 
of radio propagation, since they represent the iono- 
spheric characteristic all over the world at a fixed 
instant. The maps, given for the four hours, UT = 00, 
06, 12, 18, are presented in such a way that it is 
always midnight in LMT at the extreme left bound- 
ary and noon in LMT at the center. This type of 



presentation enables one to see the relatively small 
continuous deformations in the ionosphere as the 
earth rotates about its axis. 

Since the geographic poles are singular points in 
the numerical maps, the function r(X,^,^) is physically 
meaningless in the immediate neighborhood of the 
north and south poles. Therefore, in these regions 
dashed lines (interpolated by hand) have been in- 
serted in the contour maps to give a meaningful ex- 
ten tion to the poles. Similarly dashed lines have 
been inserted in figures 4, 5, 8, and 9. 

4.2. Polar Maps in Universal Time 

Polar maps of monthly median /o^2 iii universal 
time (Cut III), illustrated in figures 4 and 5, are 
similar to the world maps in UT (figs. 2 and 3). 
Such maps are more convenient for working with 
radio circuits passing over the polar caps than the 
(rectangular) world maps in UT shown in figures 
2 and 3. 

4.3. World Maps in Local Mean Time 

World maps of monthly median ^7^2 in local mean 
time (Cut I) represent the worldwide geographic 
variations which are not produced by changes in the 
hour angle of the sun (figs. 6 and 7). Such maps 
are particularly useful in studying the complex and, 
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SEPTEMBER, 1958 UT=00 




3o-jr-- ^-^-^-;t 



90* 120' 150* I80* 210* 240* 270* 300* 330' 360* 
= LONGITUDE 



SEPTEMBER, 1958 UT=I2 




m* 210* 240» 270* 300* 330* 0* 30* 60« 90» 120' 150* 
e = LONGITUDE 

Figure 2. World maps of monthly median foF2 (Mc/s) in 
universal time from the function T{\,6,t). 



SEPTEMBER, 1958 UT=06 



to some extent, irregular geographic variations due 
to seasonal changes, the earth's magnetic field, and 
other factors affecting the ionosphere. The maps 
are given for four different hours (LMT=00, 06, 12, 
18), illustrating the different conditions occurring at 
midnight, sunrise, noon, and late afternoon, respec- 
tively. 

4.4. Maps in Latitude and LMT 

Maps of monthly median yo^2 in latitude and LMT 
(Cut II) are given in figures 8 and 9 for several fixed 
values of longitude. These maps are similar to the 
''longitude zone" maps presently given in the 
CRPL-D Series. However, in contrast with the 
D series, it is not assumed here that longitudinal 
variation of the ionospheric characteristic is con- 
stant within a given zone. 




go* 120° 150° 



210° 240° 270° 300" 350° 0° 
9 = LONGITUDE 



Figure 3. World maps of monthly median foF2 (Mc/s) in 
universal time from the function T{\,6,t), 



4.5. Diurnal Variation 

The diurnal variation for monthly median foF^ 
(Line A) is illustrated in figures 10 and 11 for several 
different locations on the earth (the points of inter- 
section of the two meridians 90°E and 285°E with 
the five parallels 0°, ±20°, ±50°). These graphs 
illustrate the strong changes in the diurnal variation 
with geographic position. 

4.6. Longitudinal Variation (Fixed LMT) 

The longitudinal variation of monthly median 
foF2 (Line B) for a fixed instant of LMT is illustrated 
by the graphs in figure 12. The changes in this varia- 
tion (corresponding to different latitudes and hours 
of LMT) are typified by the graphs for each of the 
hours, 00 and 12, and for each of the latitudes 0°, 
±20°, ±50°. 
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NORTH POLE 
SEPTEMBER. 1958 UT=00 




SOUTH POLE 
SEPTEMBER, 1958 UT = 00 

Figure 4. Polar maps of monthly median for2 (Mc/s) in 
universal time from the function r(\,d,t). 



4.7. Latitudinal Variation (Fixed LMT) 

The latitudinal variation of monthly median foF2 
(Line C) for a fixed instant of LMT is shown by the 
graphs in figure 13. Graphs are given for each of the 
fixed longitudes 0°E, 90°E, 180°E, and 270°E, 
for each hour, LMT =00 and 12. 

4.8. Longitudinal Variation (Fixed UT) 

The variation of monthly median /oi^2 with longi- 
tude (Line D) for a fixed instant of UT is illustrated 
by the graphs in figure 14, for the two hours, UT = 00 
and 12, for each of the fixed latitudes 0°, ±20°, ±50°. 




NORTH POLE 
SEPTEMBER. 1958 UT= 12 




SOUTH POLE 
SEPTEMBER, 1958 UT=I2 

Figure 5. Polar maps of monthly median foF2 (Mc/s) in 
universal time from the function T(\,d,t). 



5. Summary 

It has been pointed out that a numerical map is 
defined by means of a relatively small table of 
numerical coefficients with which many useful appli- 
cations can be made. In particular, one can com- 
pute the value of the ionospheric characteristic at 
any desired location and instant of time. By means 
of this operation a wide variety of graphical repre- 
sentations can be computed automatically and plotted 
by machine. These include, for example, world and 
polar contour maps in either universal time or local 
mean time, illustrations of which were given in 
section 4. 
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SEPTEMBER, 1958 LMT=00 




0* 30' 



90* 120' ISO* ISO' 210' 240' 270' 300' 330' 
^ = LONGITUDE 

SEPTEMBER, 1958 LMT=I2 




0* 30' 



90* 120' 150* 180* 210' 240' 270' 300' 330 
0= LONGITUDE 



Figure 6. World yiiaps of monthly median foF2 {Mcjs) in 
local mean time from the function T{\,0,t). 



Although such graphical displays have many im- 
portant uses, it should be re-emphasized that numer- 
ical maps can be used to the greatest advantage and 
efficiency within the computer — that is, by letting 
the machine compute the values of the characteristic 
rather than by reading them manually from a contour 
map. The computer can do this much faster and 
more accurately than a person reading and inter- 
polating between maps. 

6. Appendix: Basic Problem in Computing 
Contour Maps 

Let us assume that we have been given a table of 
coefficients Dsk defining a numerical map r(X,^,/) for 
an ionospheric characteristic (for example, the 
coefficients in table 2 for monthly median /o7^2, 



SEPTEME-ZR, 1958 LMT=06 




'0* 30' 60' 90' 120' 150' 180' 210' 240 270' 300' 330* 360' 

= LONGITUDE 

SEPTEMBER, 1958 LMT= 18 




' 0* 30' 60* 



120' ISC 180* 210' 240' 270* 300* 330* W 
^= LONGITUDE 



Figure 7. World maps of monthly jnedian foF2 {Mc/s) in 
local mean time from the function V{\,d,t). 



September 1958), so that we can evaluate r(X,^,^) at 
any desired point in the three-dimensional space of 
X, 6, and t. Then the problem considered here is 
that of locating, for given fixed values of d and t, all 
values of X for which 



T(\,d,t)=lQ for some 2=0,1 



Q (16) 



(see sec. 3.2). Thus the problem consists in solving 
high degree polynomial equations in sin X for points 
X lying along selected meridians. The method we 
employ has two main steps, (1) to isolate the roots 
within certain small intervals and (2) to approach 
each root by a method of successive approximation. 
We begin the discussion with a description of two 
basic parameters which must be determined. These 
parameters are associated with the function T(\,djt), 
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SEPTEMBER, 1958 6 = 0° 




00 02 04 06 



10 12 14 16 18 20 22 24 

LOCAL MEAN TIME 



SEPTEMBER, 1958 U =180° 




00 02 04 06 08 10 12 14 16 18 20 22 24 

LOCAL MEAN TIME 

FiGUKE 8. Maps of monthly median foF2 (Mc/s) in latitude 
and LMT at fixed longitudes from the function T(\,d,t). 



6.1. Parameters e and N Associated With 
T{\d,t) 

In order to solve the problem outlined above, two 
parameters e and N are needed. The determination 
of these parameters depends upon certain general 
properties of the function T(\,d,t) as follows. The 
positive number e determines ^ the permissible error 
in an approximation to a root of (16) in the sense 
that X will be considered a solution to (16) if 



|r(XAO -7,1^6 



for some g=0,l , . . ., Q,. (17) 



For a very flat part of the function T(\,d,t), large 
variations in the independent variables make only 
small changes in T(\,d,t). Thus, for a given e, the 
error in X is proportional to the ''flatness^' of T(\,djt). 

6 As can be seen, « is net itself the upper bound of the error in X but only deter- 
mines what that upper bound will be for a particular function T{\,d,t). 



SEPTEMBER, 1958 9 = 90° 




00 02 04 06 08 10 12 14 16 

LOCAL MEAN TIME 



20 22 24 



Figure 9. Maps of monthkj median foF2 {Mc/s) in latitude 
and LMT at fixed longitudes from the function T{\,d,t). 



Since the roots of (16) are obtained by successive 
approximations, it is desirable for economy to choose 
e no larger than necessary. For most ionospheric 
characteristics a suitable choice is €=0.0001. 

The number A'^ depends upon T(\,d,t) in somewhat 
the opposite sense. As will be seen, successive 
approximations to roots of (16) are to be made by 
means of linear interpolation within the intervals 



Xn-i<X<X^ n=l,2, 



• .N 



where 



X«=[2(j,)-l]90° n=0,l,...,N. 



(18) 



(19) 



Therefore N must be chosen large enough so that, 
for fixed values of 6 and t, the function T(\,d,t) will 
be approximately linear within the intervals (18). 
If there are relative maxima or minima of T(\,d,t) 
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X= 0° 

X=-20° 

X=-50° 
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4 8 12 16 20 24 

LOCAL MEAN TIME 

Figure 10. Diurnal variation of monthly median fo¥2 from the 
function T(\,d,t). September 1958 (Fixed longitude: 
= 90°). 
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Figure 11. Diurnal variation of monthly median foF2 /rom the 
function r(X,0,O- September 1958 (Fixed longitude: 

= 285°). 
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Figure 12. Longitudinal variation of 7nonthly median foF2 
from the function T{\,6,t) for September 1958. 



within certain of the intervals (18), a few points on 
the map could be lost. 

For the sake of economy one would like to choose 
N no larger than necessary. For most ionosplieric 
characteristics it has been found that iV=36 is 
sufficient. A typical example of the variation of 
r(\,d,t) for monthly median /oi^2, September 1958, 
with 6/=270° and ^=-180° (LMT=00) is shown 
in figure 15. The vertical lines spaced every 5° 
illustrate the type of subdivision resulting from the 
choiceof A^==36. 

In the remaming part of the appendix we shall 
assume that e and A^have been appropriately chosen 
and that 6 and t are held fixed; hence for simplicity 
we shall write r(X) in place of T{\,d,t). 

6.2. Isolation of the Roots in the nth Interval 

We consider first the problem of determining the 
number of roots in the nth interval (18). We 
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Figure 13. Latitudinal variation of monthly median 10^2 
from the function V{\,d,t) for September 1958. 



define Hn as that value of g; (^=0, 1, . . ., Q) for 
which 

/..i<r(x«)^/,. (20) 

(See schematic diagram in fig. 16.) Then, since 
r (X) is approximately linear and hence monotonic 
in the nth interval, the number ^n of roots in the 
nth interval is given by 



Tn = 



\^n — ^n-l\ 



(21) 



The end points X^^ can also be roots and must 
therefore be checked. We shall denote the jn roots 
in the nth. interval by 



where 



K''' j=l, 2, 



X._i<X,^^)<X,^^)< 



Tn 



<x/^.<x. 



(22) 



(23 




LATITUDES 

50°N 

20°N 

EQUATOR 

20«»S 

SO'S 



FIXED UT=I2 



Figure 14. Longitudinal variation of monthly median foF2 ' 
from the function r{}^,e,t) for September 1958. 
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Figure 15. Latitudinal variation of monthly median foF? 
from the function V{\d,t) for September 1958. 
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and for simplicity we shall wi'ite 

r„«>=r(x„«>) i=i,2, 



T«- 



(24) 



Again, from the mono tonic behavior of r(X) in 
the nth interval it can be seen that 

(A) ju„_i</x;. impHes Tn^^^=I^^^_^+i-i 

j = l, 2, . . ., yn 

(B) Mn-i>M. implies r,^^-)=/,^_^_^ (25) 

(C) /x„_i=/i,i imphes there are no solutions in the 

71 th interval. 

Thus we have determined the number of roots in 
the Tith interval as well as the values assumed by 
r(X) at these roots. The final problem of approxi- 
mating each root successively is resolved in the 
following section. 

6.3. Method of Successive Approximation 

Since A^ was so chosen that the function r(X) is 
approximately linear within any of the intervals 
(18), a satisfactory approximation can be made to 
the roots X„/^\ X^^^^ . . ., Xn^^«^ in the ?ith interval 
by means of linear interpolation. In the procedure 
to be followed the size of the interval of interpola- 
tion is reduced after each approximation so that 
rapid convergence is assured. 

For the first root X„^^^ the initial interval of 
interpolation will be (ai^^\ bi^^^) where 



=X,^i and bi''=\ 



(26) 



so that the first approximation to X^^^^ is given by 

_6<"[r4»-r(a{»)]-a«'[r<"-r(6P')] 



rf:" 



r(6<")-r(a<») 



(27) 



(See schematic diagram in fig. 16.) In general, the 
rth approximation to tlie jtli root X„*^* is given by 



(^,"'= 



^,»)[r»>-r(a,"')]-a»'[ri^'-r(6,">)] 
r(6'«)-r(a<«) 



(28) 



where (a,'-^\ br'-^^) denotes the rth interval of interpo- 
lation for the jth root (to be defined below) . After 
each approximation we must determine which of the 
following cases occurs: 

|r<^'-r(rf<^>)|^6, (29) 

r(a<^>)<ri«<r(cZ<^-') ^ r(6<''), (30) 

r(a,'^-')gr(c?<>')<ri^'<r(6<«), (3i) 

r(a<^>)>ri^>>r«^')^ r(6;^>), (32) 

r(a<''0^r(.i;^')>ri«>r(6«). (33) 

We note that (32) and (33) cannot occur if (25A) 
holds, and (30) and (31) cannot occur if (25B) 
holds. If (29) occurs, then clearly we have the 
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Figure 10. Schematic diagram of typical nth interval. 



desired approximation, 
we let 



a4^,=ai^^ and b'/ii=d^/' 
wliereas if (31) or (:^3) occur, we let 

a'4i^d}^^ and b^r'li = K'' 



If (:]0) or (:V2) occur, then 
(34) 



(35) 



so that in any case the size of the new mterval of 
interpolation is less than that of the preceding one. 
This process is continued until (29) occurs and we 
proceed to the next root. For each root X^^^^ after 
the first (l<Ci^7n), the initial interval of interpola- 
tion is taken to be 

a[^'=K'-'' and 6p>=X„i=2, 3, . . ., 7. (36) 

so that the convergence becomes faster for each 
successive root since the initial intervals are suc- 
cessivelv reduced. 



The authors gratefully acknowledge the assistance 
they received from a number of persons. The ex- 
amples of applications of ''numerical maps,'' which 
illustrate the methods given in the present paper, 
were rendered possible by the careful work of Miss 
Martha Hinds and Miss Minadora PoKempner. 
They have been responsible for the preparation and 
analysis of the ionospheric data and for programing 
on large-scale computers. Direct assistance has 
also been received from Mrs. Estelle D. Powell, 
Mrs. Anna von Kreisler, Konald P. Graham, Mrs. 
Gladys Waggoner, and operators in the Computing 
Laboratories of the National Bureau of Standards 
(in Boulder, Colo., and in Washington, D.C.). 



661 



7. References 

Byerly, W. E., Fourier series and spherical, cylindrical, and 

ellipsoidal harmonics, pp. 144-218 (Ginn and Company, 

Boston, 1893). 
Jones, W. B., Atlas of Fourier coefficients of diurnal variation 

of foF2, NBS Tech. Note 142 (PB161643) (Apr. 1962). 
Jones, W. B., and R. M. Gallet, Ionospheric mapping by 

numerical methods, J. Intern. Telecomm. Union, No. 12, 

260-264 (Dec. 1960). 



Jones, W. B., and R. M. Gallet, The representation of diurnal 
and geographic variations of ionospheric data by numerical 
methods. Telecommunication J. 29, No. 5, 129-149 (May 
1962), and J. Research NBS 66D (Radio Prop.) No. 4, 
419-438 (July- Aug. 1962). 

Wright, J. W., L. R. Wescott, and D. J. Brown, Mean electron 
density variations of the quiet ionosphere 4 — June 1959, 
NBS Tech. Note 49-4 (May 1961). 

(Paper 66D6-225) 



662 



